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Double beta decay is indispensable to solve the question of the neutrino mass 
matrix together with v oscillation experiments. Recent analysis of the most sen- 
sitive experiment since nine years - the HEIDELBERG-MOSCOW experiment in 
Gran-Sasso - yields a first indication for the neutrinolcss decay mode. This result 
is the first evidence for lepton number violation and proves the neutrino to be 
a Majorana particle. We give the present status of the analysis in this report. 
It excludes several of the neutrino mass scenarios allowed from present neutrino 
oscillation experiments - only degenerate scenarios and those with inverse mass 
hierarchy survive. This result allows neutrinos to still play an important role as 
dark matter in the Universe. To improve the accuracy of the present result, con- 
siderably enlarged experiments are required, such as GENIUS. A GENIUS Test 
Facility has been funded and will come into operation by early 2003. 



1. Introduction 

Double beta decay is the most sensitive probe to test lepton number con- 
servation. Further it seems to be the only way to decide about the Dirac or 
Majorana nature of the neutrino. Observation of Of/3/3 decay would prove that 
the neutrino is a Majorana particle and would be another clear sign of beyond 
standard model physics. Recently atmospheric and solar neutrino oscillation 
experiments have shown that neutrinos are massive. This was the first indi- 
cation of beyond standard model physics. The absolute neutrino mass scale, 
however, cannot be determined from oscillation experiments alone. Double 
beta decay is indispensable also to solve this problem. 

The observable of double beta decay is the effective neutrino mass 
(to) = | J2 U^rrii | = \m£ | + |ml? | + e^ 3 \m$ | , 
with U ei denoting elements of the neutrino mixing matrix, rrii neutrino mass 
eigenstates, and 4>i relative Majorana CP phases. It can be written in terms 
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of oscillation parameters 5 



|m£>| = lUe^m, 

\m^\ = \U e2 \ 2 ^Am 2 21 +m 2 , 



(2) 



(1) 



m^\ = \U e3 \ 2 ^Am 2 32 + Am 2 21 +ml 



(3) 



The effective mass (m) is related with the half-life for 0v(3(3 decay via 

(■^1/2) ~ ( m ^) 2 i an( l f° r the limit on Ty 2 dcduciblc in an experiment we 
have 



Here a is the isotopical abundance of the (3(3 emitter; M is the active 
detector mass; t is the measuring time; AE is the energy resolution; B is 
the background count rate and e is the efficiency for detecting a (3(3 signal. 
Determination of the effective mass fixes the absolute scale of the neutrino 
mass spectrum. 5,T ). 

In this paper we will discuss the status of double beta decay search. We 
shall, in section 2, discuss the recent evidence for the neutrinoless decay mode, 
from the HEIDELBERG-MOSCOW experiment and the consequences for the 
neutrino mass scenarios which could be realized in nature. In section 3 we 
discuss the possible future potential of 0v(3(3 experiments, which could improve 
the present accuracy. 

2. Evidence for the Neutrinoless Decay Mode 

The status of present double beta experiments is shown in Fig. 1 and is ex- 
tensively discussed in 7 . The HEIDELBERG-MOSCOW experiment using the 
largest source strength of 11 kg of enriched 76 Ge (enrichment 86%) in form of 
five HP Ge-detectors is running since August 1990 in the Gran-Sasso under- 
ground laboratory 7 > 3 > 4 , 1,32,28 ^ anc j j g smce nme years now the most sensitive 
double beta experiment worldwide. 

2.1. Data from the HEIDELBERG-MOSCOW Experiment 

The data taken in the period August 1990 - May 2000 (54.9813 kg y, or 
723.44 mol-years) are shown in Fig. 2 in the section around the Qpp value 
of 2039.006 keV 16 ' 17 '. Fig. 2 is identical with Fig. 1 in \ except that we show 
here the original energy binning of the data of 0.36 keV. These data have been 
analysed 1,2,3 with various statistical methods, with the Maximum Likelihood 
Method and in particular also with the Bayesian method (see, e.g. n ). This 
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Figure 1. Present sensitivity, and expectation for the future, of the most promising j3/3 
experiments. Given are limits for (m), except for the HEIDELBERG-MOSCOW experiment 
where the recently observed value is given (95% c.l. range and best value). Framed parts of 
the bars: present status; not framed parts: future expectation for running experiments; solid 

and dashed lines: experiments under construction or proposed, respectively. For references 

7,2,3,35,46 



method is particularly suited for low counting rates, where the data follow a 
Poisson distribution, that cannot be approximated by a Gaussian. Details and 
the results of the analysis are given in 1 < 2 > 3 . 

Our peak search procedure (for details see 2 ' 3 ) reproduces (see li2 ' 3 ) by- 
lines at the positions of known weak lines from the decay of 214 Bi at 2010.7, 
2016.7, 2021.8 and 2052.9 keV 15 . In addition, a line centered at 2039 keV 
shows up (see Fig. 3). This is compatible with the Q-value 16,17 of the double 
beta decay process. The Bayesian analysis yields, when analysing a ±5<r range 
around Qpp (which is the usual procedure when searching for resonances in 
high-energy physics) a confidence level (i.e. the probability K) for a line to 
exist at 2039.0 keV of 96.5 % c.l. (2.1 a) (see Fig. 3). We repeated the 
analysis for the same data, but except detector 4, which had no muon shield 
and a slightly worse energy resolution (46.502 kg y). The probability we find 
for a line at 2039.0 keV in this case is 97.4% (2.2 a) 1 ' 2 ' 3 . 

Fitting a wide range of the spectrum yields a line at 2039 keV at 91% c.l. 
(see Fig. 2). 

We also applied the Feldman-Cousins method 12 . This method (which does 
not use the information that the line is Gaussian) finds a line at 2039 keV on 
a confidence level of 3.1 a (99.8% c.l.). In addition to the line at 2039 keV 
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Figure 2. The spectrum taken with the 76 Ge detectors Nr. 1,2,3,4,5 over the period August 
1990 - May 2000 (54.9813 kg y) in the original 0.36 keV binning, in the energy range 2000 - 
2100 keV. Simultaneous fit of the 214 Si lines and the two high-energy lines yield a probability 
for a line at 2039.0 keV of 91%. 
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Figure 3. Left: Probability K that a line exists at a given energy in the range of 2000-2080 
keV derived via Bayesian inference from the spectrum shown in Fig. 2. Right: Result of 
a Bayesian scan for lines as in the left part of this figure, but in an energy range of ±5cr 
around Q33. 

we find candidates for lines at energies beyond 2060 keV and around 2030 keV, 
which at present cannot be attributed. This is a task of nuclear spectroscopy. 

Important further information can be obtained from the time structures of 
the individual events. Double beta events should behave as single site events 
(see Fig. 4 left), i.e. clearly different from a multiple scattered 7-event (see Fig. 
4 right). It is possible to differentiate between these different types of events 
by pulse shape analysis. We have developped three methods of pulse shape 
analysis 8 ' 9 > 10 during the last seven years, one of which has been patented and 
therefore only published recently. 

Installation of Pulse Shape Analysis (PSA) has been performed in 1995 for 
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Figure 4. Left: Shape of one candidate for 0vf3f3 decay classified as SSE by all three 
methods of pulse shape discrimination. Right: Shape of one candidate classified as MSE by 
all three methods. 

the four large detectors. Detector Nr. 5 runs since February 1995, detectors 
2,3,4 since November 1995 with PSA. The measuring time with PSA from 
November 1995 until May 2000 is 36.532 kg years, for detectors 2,3,5 it is 
28.053 kg y. 

In the SSE spectrum obtained under the restriction that the signal simul- 
taneously fulfills the criteria of all three methods for a single site event, we find 
again indication of a line at 2039.0 keV (see Figs. 5, 6). 

We find 9 SSE events in the region 2034.1 - 2044.9keV (± 3er around Q^). 
Bayes analysis of the range 2032 - 2046 keV yields a signal of single site events, 
as expected for neutrinoless double beta decay, with 96.8% c.l. at the 
value. The Feldman-Cousins method gives a signal at 2039.0 keV of 2.8 a 
(99.4%). 

The analysis of the line at 2039.0 keV before correction for the efficiency 
yields 4.6 events (best value) or (0.3 - 8.0) events within 95% c.l. ((2.1 - 
6.8) events within 68.3% c.l.). Corrected for the efficiency to identify an SSE 
signal by successive application of all three PSA methods, which is 0.55 ± 0.10, 
we obtain a Qvj3(3 signal with 92.4% c.l.. The signal is (3.6 - 12.5) events with 
68.3% c.l. (best value 8.3 events). Thus, with proper normalization concerning 
the running times (kgy) of the full and the SSE spectra, we see that almost 
the full signal remains after the single site cut (best value), while the 214 Bi 
lines (best values) are considerably reduced. We have used a 238 Th source to 
test the PSA method. We find the reduction of the 2103 keV and 2614 keV 
228 Th lines (known to be multiple site or mainly multiple site), relative to the 
1592 keV 22S Th line (known to be single site), shown in Fig. 7. This proves 
that the PSA method works efficiently. Essentially the same reduction as for 
the Th lines at 2103 and 2614 keV and for the weak Bi lines is found for the 
strong 214 Bi lines (e.g. at 609.6 and 1763.9keV (Fig. 7)). 

The possibility, that the single site signal is the double escape line cor- 
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Figure 5. Scan for lines in the single site event spectrum taken from 1995-2000 with de- 
tectors Nr. 2,3,5, (Fig. 6), with the Bayesian method. Left: Energy range 2000 -2080 keV. 
Right: Energy range of analysis ±4.4cr around Q^. 
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Figure 6. Sum spectrum of single site events, measured with the detectors Nr. 2,3,5 op- 
erated with pulse shape analysis in the period November 1995 to May 2000 (28.053 kg y), 
summed to IkeV bins. Only events identified as single site events (SSE) by all three pulse 
shape analysis methods 8 " 9 > 10 have been accepted. The curve results from Bayesian inference 
in the way explained in sec. 3. When corrected for the efficiency of SSE identification (see 
text), this leads to the following value for the half-life: T°^ 2 =(0.88 - 22.38)x 10 25 y (90% 
c.l.). 

responding to a (much more intense!) full energy peak of a 7-line, at 
2039+1022=3061keV is excluded from the high-energy part of our spectrum. 

2.2. Comparison with earlier results 

We applied the same methods of peak search as used in our analysis to the 
spectrum, measured in the Ge experiment by Caldwell et al. 18 more than a 
decade ago. These authors had the most sensitive experiment using natural Ge 
detectors (7.8% abundance of 76 Ge). With their background being a factor of 
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Figure 7. Relative suppression ratios: Remaining intensity after pulse shape analysis com- 
pared to the intensity in the full spectrum. Right: Result of a calibration measurement with 
a Th source - ratio of the intensities of the 1592 keV line (double escape peak, known to 
be 100% SSE), set to 1. The intensities of the 2203 keV line (single escape peak, known to 
be 100% MSE) are strongly reduced (error bars are ±lo\ The same order of reduction is 
found for the strong Bi lines occuring in our spectrum - shown in this figure are the lines 
at 609.4 and 1763.9 keV. Left: The lines in the range of weak statistics around the line at 
2039 keV (shown are ratios of best fit values). The Bi lines are reduced compared to the line 
at 2039 keV (set to 1), as to the 1592 keV SSE Th line. 

9 higher than in the present experiment, and their measuring time of 22.6 kg y, 
they have a statistics for the background larger by a factor of almost 4 in their 
(very similar) experiment. This allows helpful conclusions about the nature of 
the background. 

The peak scanning finds 3 (Fig. 8) indications for peaks essentially at the 
same energies as in Fig. 3. This shows that these peaks are not fluctuations. 
In particular it sees the 2010.78, 2016.7, 2021.6 and 2052.94keV ni Bi lines. 
It finds, however, no line at Qpp. This is consistent with the expectation 
from the rate found from the HEIDELBERG-MOSCOW experiment. About 
16 observed events in the latter correspond to 0.6 expected events in the Cald- 
well experiment, because of the use of non-enriched material and the shorter 
measuring time. 

Another Ge experiment (IGEX) using 8.8 kg of enriched 76 Ge, but collect- 
ing since beginning of the experiment in the early nineties till shutdown in end 
of 1999 only 8.8 kg y of statistics 42 , because of this low statistics also natu- 
rally cannot see any signal at 2039 keV. This is consistent with our result. The 
number of counts expected in that experiment from our result would be about 
only 2.6. In their spectrum however, which is unfortunately shown only in the 
energy range 2020-2060 keV, the line around 2030 keV, indicated in Figs. 3,6, 
is also clearly seen. 
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Figure 8. Peak scanning of the spectrum measured by Caldwell et al. 8 , with the Maximum 
Likelihood method (upper part), and with the Baycsian method (lower part) (as in Figs. 3,6) 
(see 3 ). 



2.3. Proofs and Disproofs 

The result described in section 2.1 has been questioned in some papers (Aalseth 
et al, hep-ex/0202018; Feruglio et al., Nucl. Phys. B 637(2002)345; Zdesenko 
et al., Phys. Lett. B 546(2002) 206, and Kirpichnikov, talk at Meeting of 
Physical Section of Russian Academy of Sciences, Moscow, December 2, 2002, 
and priv. communication, Dec. 3, 2002.) These claims against our results are 
incorrect in various ways. 

The arguments in the first two of these papers can be easily rejected. We 
have published this in hep-ph/0205228. In particular their estimates of the in- 
tensities of the 214 Bi lines are simply wrong, because the summing effect of the 
energies of consecutive gamma lines (more precisely called True Coincidence 
Summing - TCS) which is an elementary effect known in nuclear spectroscopy, 
was not taken into account (for a detailed discussion see 3 ). Further none of 
the papers has performed a Monte Carlo simulation of our setup, which is the 
only way to come to quantitative statements. 

The paper by Zdesenko et al. starts from an arbitrary assumption, namely 
that there are lines in the spectrum at best only at 2010 and 2053 keV. This 
contradicts to the experimental result, according to which there are further 
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lines in the spectrum (see Figs. 3,5,8 in the present paper). In this way 
they artificially increase the background in their analysis and come to wrong 
conclusions. 

All three of these papers, when discussing the choice of the width of the 
search window, ignore the results of the statistical simulations we give in hep- 
ph/0205228 and we have published in 2 and 3 . 

Kirpichnikov states that from his analysis he clearly sees the 2039 keV line 
in the full (not pulse shape discriminated) spectrum with > 99% c.l. He claims 
that he does not see the signal in the pulse shape spectrum. One reason to 
see less intensity certainly is that in this case he averages for determination of 
the background over the full energy range without allowing for any lines. His 
result is in contradiction with the result we obtain under the same assumption 
of assuming just one line (at Q/3,3) and a continuous background (see Fig. 5 of 
this paper). 

2.4. Half-Life and Effective Neutrino Mass 

Having shown that the signal at consists of single site events and is not 

a 7-line, we translate the observed number of events into half-lifes. We obtain 
1,2,3 T o^ = ( 8 _ lg 3) x 1Q 25 y ( 95% c j ) with a best value of L5 x 1Q 25 y 

Assuming that the §v(3[3 amplitude is dominated by the neutrino mass mech- 
anism, we obtain, with the nuclear matrix element from 13 an effective mass 
of {m v )=(0.ll - 0.56) cV (95% c.l.). 

The result obtained is consistent with all other double beta experiments - 
which still reach less sensitivity. The most sensitive experiments following the 
HEIDELBERG-MOSCOW experiment arc the geochemical 128 Te experiment 
with T^ 2 > 2(7.7) x 10 24 y (68% c.l.), 20 the 136 le experiment by the DAMA 
group with T\ v /2 > 1.2 x 10 24 y (90% c.l.) 21 , a second 76 Ge experiment with 
T°J 2 > 1.2 x 10 24 y 19 and a nat Ge experiment with T°^ 2 > 1 x 10 24 y 18 < 19 . 
Other experiments are already about a factor of 100 less sensitive concerning 
the 0v/3f3 half-life: the Gotthard TPC experiment with 136 Ae yields 22 T°^ 2 > 
4.4 x 10 23 y (90% c.l.) and the Milano Mibeta cryodetector experiment 44 
T^ 2 > 1.44 x 10 23 y (90% c.l.). 

Another experiment 42 with enriched 76 Ge, which has stopped operation in 
1999 after reaching a significance of 8.8 kg y, yields (if one believes their method 
of 'visual inspection' in their data analysis), in a conservative analysis, a limit 
of about T°J 2 > 5 x 10 24 y (90% c.l.). The 128 Te geochemical experiment 
yields (m„) < 1.1 eV (68 % c.l.) 20 , the DAMA 136 Xe experiment (m„) < 
(1.1 - 2.9) eV 21 and the 130 Te cryogenic experiment yields (m v ) < 1.8 eV 44 . 

Concluding we obtain, with about 95% probability, first evidence for the 
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neutrinoless double beta decay mode. As a consequence, at this confidence 
level, lepton number is not conserved. Further the neutrino is a Majorana 
particle. If the Qv(3(3 amplitude is dominated by exchange of a massive neutrino 
the effective mass (m) is deduced to be (m) = (0.11 - 0.56) eV (95% c.L), with 
best value of 0.39 eV. Allowing conservatively for an uncertainty of the nuclear 
matrix elements of ± 50% (for detailed discussions of the status of nuclear 
matrix elements we refer to 7 > 3 ' 14 and references therein) this range may widen 
to (m) = (0.05 - 0.84) eV (95% c.L). 

Assuming other mechanisms to dominate the Qi/(3(3 decay amplitude, the 
result allows to set stringent limits on parameters of SUSY models, lepto- 
quarks, compositcncss, masses of heavy neutrinos, the right-handed W boson 
and possible violation of Lorentz invariance and equivalence principle in the 
neutrino sector. For a discussion and for references we refer to 7 > 2 3,26,34,46_ 

With the limit deduced for the effective neutrino mass, the HEIDELBERG- 
MOSCOW experiment excludes several of the neutrino mass scenarios allowed 
from present neutrino oscillation experiments - allowing only for degenerate 
and inverse hierarchy mass scenarios 6 . 

Assuming the degenerate scenarios to be realized in nature we fix - ac- 
cording to the formulae derived in 5 - the common mass eigenvalue of the 
degenerate neutrinos to m = (0.05 - 3.4) eV. Part of the upper range is already 
excluded by tritium experiments, which give a limit of m < 2.2-2.8 eV (95% 
c.L) 30 . The full range can only partly (down to ~ 0.5 eV) be checked by 
future tritium decay experiments, but could be checked by some future (3(3 
experiments (see, e.g., next section). The deduced best value for the mass is 
consistent with expectations from experimental fi — > 67 branching limits in 
models assuming the generating mechanism for the neutrino mass to be also 
responsible for the recent indication for as anomalous magnetic moment of 
the muon 50 . It lies in a range of interest also for Z-burst models recently 
discussed as explanation for super-high energy cosmic ray events beyond the 
GKZ-cutoff 48 > 49 ' 52 . A recent model with underlying A4 symmetry for the 
neutrino mixing matrix also leads to degenerate neutrino masses consistent 
with the present result from 0v(3(3 decay 51 . The range of (m) fixed in this 
work is, already now, in the range to be explored by the satellite experiments 
MAP and PLANCK 27 - 5 (see Fig. 10). 

The neutrino mass deduced leads to 0.002> fl„h 2 < 0.1 and thus may allow 
neutrinos to still play an important role as hot dark matter in the Universe 33 . 

3. Future of (3(3 Experiments - GENIUS and Other Proposals 

With the HEIDELBERG-MOSCOW experiment, the era of the small smart 
experiments is over. New approaches and considerably enlarged experiments 
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Figure 9. The impact of the evidence obtained for neutrinoless double beta decay in this 
paper (best value of the effective neutrino mass (m) = 0.39 eV, 95% confidence range (0.05 
- 0.84) eV - allowing already for an uncertainty of the nuclear matrix element of a factor 
of ± 50%) on possible neutrino mass schemes. The bars denote allowed ranges of (m) in 
different neutrino mass scenarios, still allowed by neutrino oscillation experiments (see 6 ). 
Hierarchical models are excluded by the new 0vf3f3 decay result. Also shown are the expected 
sensitivities for the future potential double beta experiments CUORE, MOON, EXO and 
the 1 ton and 10 ton project of GENIUS 7,34,47,25 



(as discussed, e.g. in 7 ' 26, 47, 25 ' 29 ) will be required in future to fix the neutrino 
mass with higher accuracy. 

Since it was realized in the HEIDELBERG-MOSCOW experiment, that the 
remaining small background is coming from the material close to the detector 
(holder, copper cap, ...), elimination of any material close to the detector 
will be decisive. Experiments which do not take this into account, like, e.g. 
CUORE 43 , and MAJORANA 45 , will allow at best only rather limited steps in 
sensitivity. Furthermore there is the problem in cryodetectors that they cannot 
differentiate between a (3 and a 7 signal, as this is possible in Ge experiments. 

Another crucial point is - see eq. (4) - the energy resolution, which can 
be optimized only in experiments using Germanium detectors or bolometers. 
It will be difficult to probe evidence for this rare decay mode in experiments, 
which have to work - as result of their limited resolution - with energy win- 
dows around of several hundreds of keV, such as NEMO III 41 , EXO 39 , 
CAMEO 40 . 

Another important point is (see eq. 4), the efficiency of a detector for detec- 



12 



mi (eV) HEIDELBERG-MOSCOW Positive EVIDENCE 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



tan 2 6 

Figure 10. Double beta decay observable (m) and oscillation parameters: The case for 
degenerate neutrinos. Plotted on the axes are the overall scale of neutrino masses mo and 
mixing tan 2 9 12 . Also shown is a cosmological bound deduced from a fit of CMB and 
large scale structure 27 and the expected sensitivity of the satellite experiments MAP and 
PLANCK. The present limit from tritium f3 decay of 2.2 eV 30 would lie near the top of the 
figure. The range of (m) fixed by the HEIDELBERG-MOSCOW experiment 1 is, in the case 
of small solar neutrino mixing, already in the range to be explored by MAP and PLANCK 
27 



tion of a f3(3 signal. For example, with 14% efficiency a potential future 100 kg 
82 Se NEMO experiment would be, because of its low efficiency, equivalent only 
to a 10kg experiment (not talking about the energy resolution). 

In the first proposal for a third generation double beta experiment, the 
GENIUS proposal ^i 24 * 26 * 47 ,^ the idea is to use 'naked' Germanium detec- 
tors in a huge tank of liquid nitrogen. It seems to be at present the only 
proposal, which can fulfill both requirements mentioned above - to increase the 
detector mass and simultaneously reduce the background drastically. GENIUS 
would - with only 100 kg of enriched 76 67e - increase the confidence level of the 
present pulse shape discriminated Ov/3/3 signal to 4cr within one year, and to 7a 
within three years of measurement (a confirmation on a 4cr level by the MA- 
JORANA project would at least need ^230 years, the CUORE project would 
need (ignoring for the moment the problem of identification of the signal as 
a p(3 signal) 3700 years). With ten tons of enriched 76 Ge GENIUS should 
be capable to investigate also whether the neutrino mass mechanism or an- 
other mechanism (see, e.g. 7 ) is dominating the Qufifi decay amplitude. A 
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GENIUS Test Facility is at present under construction in the GRAN SASSO 
Underground Laboratory 37 > 36 . 

4. Conclusion 

The status of present double beta decay search has been discussed, and re- 
cent evidence for a non-vanishing Majorana neutrino mass obtained by the 
HEIDELBERG-MOSCOW experiment has been presented. Future projects to 
improve the present accuracy of the effective neutrino mass have been briefly 
discussed. The most sensitive of them and perhaps at the same time most 
realistic one, is the GENIUS project. GENIUS is the only of the new projects 
which simultaneously has a huge potential for cold dark matter search, and for 
real-time detection of low-energy neutrinos (see 23 > 29 > 31 > 34 > 38 . 35 > 46 ) 

5. Acknowledgment 

The author is indebted to his colleaques A. Dietz and I.V. Krivosheina for 
the fruitful and pleasant collaboration, and to C. Tomei and C. Dorr for their 
help in the analysis of the Bi lines. The generous hospitality by Profs. D. 
Ahluwalia and M. Kirchbach during the Zacatecas Forum on Physics is grate- 
fully acknowleged. 

References 

1. H.V. Klapdor-Kleingrothaus et al. hep-ph/0201231 and Mod. Phys. Lett. A 16 

(2001) 2409-2420; 

2. H.V. Klapdor-Kleingrothaus, A. Dietz and I.V. Krivosheina, Part, and Nucl. 110 

(2002) 57-79. 

3. H.V. Klapdor-Kleingrothaus, A. Dietz and I.V. Krivosheina, Foundations of 
Physics 31 (2002) 1181-1223. 

4. H.V. Klapdor-Kleingrothaus et al., (HEIDELBERG-MOSCOW Collaboration), 
Eur. Phys. J. A 12 (2001) 147 and hep-ph/0103062, and in Proc. of "Third 
Int. Conf. on Dark Matter in Astro- and Part. Phys.", DARK2000 H.V. 
Klapdor-Kleingrothaus (Editor), Springer- Verlag Heidelberg, (2001) 520 - 533 
and http : / /www.mpi — hd.mpg.de/non_acc/. 

5. H.V. Klapdor-Kleingrothaus, H. Pas and A.Yu. Smirnov, Phys. Rev. D 63 (2001) 
073005 and hep-ph/0003219; in Proc. of DARK'2000, Heidelberg, 10-15 July, 
2000, Germany, ed. H.V. Klapdor-Kleingrothaus, Springer, Heidelberg (2001) 420. 

6. H.V. Klapdor-Kleingrothaus and U. Sarkar, Mod. Phys. Lett. A 16 (2001) 2409 

7. H.V. Klapdor-Kleingrothaus, "60 Years of Double Beta Decay - From Nuclear 
Physics to Beyond the Standard Model", World Scientific, Singapore (2001) 1281 p. 

8. J. Hellmig and H.V. Klapdor-Kleingrothaus, Nucl. Instrum. Meth. A 455 (2000) 
638 - 644. 

9. J. Hellmig, F. Petry and H.V. Klapdor-Kleingrothaus, Patent DE19721323A. 

10. B. Majorovits and H.V. Klapdor-Kleingrothaus. Eur. Phys. J. A 6 (1999) 463. 



14 



11. G. D'Agostini, hep-ex/0002055, W. von der Linden and V. Dose, Phys. Rev. 
E 59 6527 (1999), and F.H. Frohner, JEFF Report 18 NEA OECD (2000) 
and Nucl. Sci. a. Engineering 126 (1997) 1, K. Weise and W. Woger, "Mes- 
sunsicherheit und Messdatenauswertung" , Wiley- VCH, Weinheim 1999, and P.M. 
Lee, "Bayesian Statistics: An Introduction", Second edition, Arnold, London 1997, 
and A. O'Hagan, "Bayesian Inference" , Volume 2B of "Kendall's Advanced Theory 
of Statistics", Arnold, London, 1997. 

12. D.E Groom ct al., Particle Data Group, Eur. Phys. J. C 15 (2000) 1; G.J. 
Feldman and R.D. Cousins, Phys. Rev. D 57 (1998) 3873. 

13. A. Staudt, K. Muto and H.V. Klapdor-Kleingrothaus, Eur. Lett. 13 (1990) 31. 

14. T. Tomoda, Kept. Prog. Phys. 54 (1991) 53 - 126. 

15. R.B. Firestone and V.S. Shirley, Table of Isotopes, 8-th Edition, John Wiley and 
Sons, Incorp., N.Y. (1998). 

16. G. Douysset et al., Phys. Rev. Lett. 86 (2001) 4259 - 4262. 

17. J.G. Hykawy et al., Phys. Rev. Lett. 67 (1991) 1708. 

18. D. Caldwell, J. Phys. G 17, S137-S144 (1991). 

19. A. A. Vasenko, I.V. Kirpichnikov et al., Mod. Phys. Lett.A 5 (1990) 1299-1306. 

20. O. Manuel et al., in Proc. Intern. Conf. Nuclear Beta Decays and the Neutrino, 
eds. T. Kotani et al., World Scientific (1986) 71, J. Phys. G: Nucl. Part. Phys. 
17 (1991) S221-S229; T. Bernatovicz et al. Phys. Rev. Lett. (1992) 2341. 

21. R. Bernabei et al., Phys. Lett. (2002) 23-28. 

22. R. Luscher et al., Phys. Lett. (1998) 407. 

23. H.V. Klapdor-Kleingrothaus in Proceedings of BEYOND'97, First International 
Conference on Particle Physics Beyond the Standard Model, Castle Ringberg, 
Germany, 8-14 June 1997, edited by H.V. Klapdor-Kleingrothaus and H. Pas, 
IOP Bristol (1998) 485 - 531 and in Int. J. Mod. Phys. A 13 (1998) 3953. 

24. H.V. Klapdor-Kleingrothaus, J. Hellmig and M. Hirsch, J. Phys. G 24 (1998) 
483-516. 

25. H.V. Klapdor-Kleingrothaus et al. MPI-Report MPI-H-V26-1999 and 
Preprint: hep-ph/9910205 and in Proceedings of the 2nd Int. Conf. on Particle 
Physics Beyond the Standard Model BEYOND'99, Castle Ringberg, Germany, 6- 
12 June 1999, edited by H.V. Klapdor-Kleingrothaus and I.V. Krivosheina, IOP 
Bristol (2000) 915 - 1014. 

26. H.V. Klapdor-Kleingrothaus, in Proc. of 18th Int. Conf. on NEUTRINO 98, 
Takayama, Japan, 4-9 Jun 1998, (eds) Y. Suzuki et al. Nucl. Phys. Proc. Suppl. 
77 (1999) 357 - 368. 

27. R.E. Lopez, astro-ph/9909414; J.R. Primack and M.A.K. Gross, astro- 
p/i/ 0007165 ; J- r - Primack, astro-ph/0007187; J. Einasto, in Proc. of 
DARK2000, Heidelberg, Germany, July 10-15, 2000, Ed. H.V. Klapdor- 
Kleingrothaus, Springer, Heidelberg, (2001). 

28. HEIDELBERG-MOSCOW Coll. (M. Gunther et al.), Phys. Rev. D55(1997)54. 

29. H.V. Klapdor-Kleingrothaus, in Proc. of Int. Conference NOW2000 - "Origins 
of Neutrino Oscillations", ed. G. Fogli, Nucl. Phys. B 100 (2001) 309 - 313. 

30. J. Bonn et al., in Proc. of NEUTRINO2000, (2000), ed. J. Law et al. Nucl. Phys. 
B 91 (2001) 273 - 279. 

31. V.A. Bednyakov and H.V. Klapdor-Kleingrothaus, /iep-p/i/0011233 and Phys. 
Rev. D 63 (2001) 095005. 



15 



32. H.V. Klapdor-Kleingrothaus, in Proc. of the Int. Symposium on Advances in 
Nuclear Physics, eds.: D. Poenaru and S. Stoica, World Scientific, Singapore 
(2000) 123 - 129. 

33. H.V. Klapdor-Kleingrothaus, Int. J. Mod. Phys. A17 (2002) 3421-3431, in Proc. 
of Intern Conf. LP01, Rome, Italy, July 2001. 

34. H.V. Klapdor-Kleingrothaus, Springer Tracts in Modern Physics, 163 (2000) 
69-104, Springer-Verlag, Heidelberg, Germany (2000). 

35. H.V. Klapdor-Kleingrothaus, in Proc. Int. Workshop on Low Energy Solar Neu- 
trinos, LowNu2, Dec. 4-5. ed: Y. Suzuki et al. World Scientific, Singapore (2001). 

36. H.V. Klapdor-Kleingrothaus, Internal Report MPI-H-V32-2000. 

37. H.V. Klapdor-Kleingrothaus et al., hep-ph/0103082 , Nucl. Instr. Meth. A 481 
(2002) 149-159. 

38. H.V. Klapdor-Kleingrothaus and I.V. Krivosheina, in these Proc. 

39. G. Gratta, talk given on ApPEC, Astroparticle Physics European Coordination, 
Paris, France 22.01.2002, and in Proc. of LowNu2, Dec. 4-5 (2000) Tokyo, Japan, 
ed: Y. Suzuki, World Scientific (2001) p.98. 

40. CAMEO Collaboration in Proc. of Taup 2001, Gran Sasso, Italy, 8.-12. Sept. 
2001 (editors A. Bettini et al.). 

41. X. Sarazin (for the NEMO Collaboration), talk given on ApPEC, Astroparticle 
Physics European Coordination, Paris, France 22.01.2002; (NEMO Collabora- 
tion), Contr. paper for XIX Int. Conf. NEUTRINO2000, Sudbury, Canada, June 
16 - 21, 2000 LAL 00-31 (2000) 1 - 10 and (NEMO-III Collaboration) in Proc. 
of NANPino2000, Dubna, Russia, July 2000, ed. V. Bednjakov, Part, and Nucl. 
Lett. 3, 62 (2001). 

42. C.E. Aalseth et al. (IGEX Collaboration), Yad. Fiz. 63, No 7, 1299 - 1302 
(2000); Phys. Rev. D 65 (2002) 092007. 

43. A. Giuliani et al., Nucl. Phys. Proc. Suppl. 110 (2002) 64-66. 

44. A. Alessandrello et al., Phys. Lett. B 486 (2000) 13-21 and S. Pirro et al., Nicl. 
Instr. Methods A 444 (2000) 71-76. 

45. L. DeBraekeleer, talk at Worksh. on the Next Generation U.S. Underground 
Science Facility, WIPP, June 12-14, 2000, Carlsbad, New Mexico, USA; C.E. 
Aalseth et al. (Majorana Collaboration), in Proc. of TAUP'2001, Gran Sasso, 
Italy, September 2001, ed. A. Bettini et al. (2002). 

46. H.V. Klapdor-Kleingrothaus, in Proc. of NANPino-2000, Dubna, July 19-22, 
2000, Particles and Nuclei, Letters iss. 1/2 (2001) and hep-ph/ 0102319. 

47. H.V. Klapdor-Kleingrothaus, hep-ph/0103074 and in Proc. of NOON 2000, 
Tokyo, Japan, 6-18 Dec. 2000, ed: Y. Suzuki et al., World Scientific (2001). 

48. D. Fargion et al., in Proc. of DARK2000, Heidelberg, Germany, July 10-15, 2000, 
Ed. H.V. Klapdor-Kleingrothaus, Springer, Heidelberg, (2001) 455 - 468, and in 
Proc. of Beyond the Desert 2002, BEYOND02, Oulu, Finland, June 2002, IOP 
2003, ed. H.V. Klapdor-Kleingrothaus. 

49. T.J. Weiler, in Proc. Beyond the Desert 1999, Tegernsee, Germany, 6-12 Juni 
1999, edited by H.V. Klapdor-Kleingrothaus and I.V. Krivosheina, IOP Bristol 

(2000) 1085 - 1106; H. Pas and T.J. Weiler, Phys. Rev. D 63 (2001) 113015. 

50. E. Ma and M. Raidal, Phys. Rev. Lett. 87 (2001) 011802; Erratum-ibid. 87 

(2001) 159901. 

51. K.S. Babu, E. Ma and J.W.F. Valle, hep-ph/0206292. 



16 



52. Z. Fodor et al., JHEP (2002) 0206:046, or hep-ph/0203198, and Z. Fodor, in Proc. 
of Intern. Conf. on Physics Beyond the Standard Model: Beyond the Desert 02, 
BEYOND'02, Oulu, Finland, 2-7 Jun 2002, IOP, Bristol, 2003, ed. H.V. Klapdor- 
Kleingrothaus and hep-ph/0210123. 



